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ABSTRACT 

Researchers at Los Alamos National Laboratory are 
considering the possibility of using the Low Energy 
Demonstration Accelerator, constructed at Los Alamos 
Neutron Science Center for the Accelerator Production 
of Tritium Project, as a neutron source. Evaluated nu- 
clear data are needed for the p+^Li reaction, to predict 
neutron production from thin and thick lithium tar- 
gets. In this paper we describe evaluation methods that 
make use of experimental data, and nuclear model cal- 
culations, to develop an ENDF- formatted data library 
for incident protons with energies up to 150 MeV. The 
important ^Li(p,no) and ^Li(p,ni) reactions are eval- 
uated from the experimental data, with their angular 
distributions represented using Lengendre polynomial 
expansions. The decay of the remaining reaction flux 
is estimated from GNASH nuclear model calculations. 
This leads to the emission of lower-energy neutrons and 
other charged particles and gamma-rays from preequi- 
librium and compound nucleus decay processes. 

The evaluated ENDF-data are described in detail, 
and illustrated in numerous figures. We also illustrate 
the use of these data in a representative application by a 
radiation transport simulation with the code MCNPX. 

L INTRODUCTION 

As a part of the Accelerator Production of Tritium 
(APT) Project [1], the Low Energy Demonstration 
Accelerator (LEDA) has been constructed at the Los 
Alamos Neutron Science Center (LANSCE) [2]. LEDA 
is a high-current low-energy proton accelerator that can 
be used to provide a source of neutrons, following pro- 
ton bombardment on suitable targets. For instance, 
high-Z targets can be used to produce spallation neu- 
trons. However, there is a recent interest in the use 
of a ^Li target, which, when bombarded with protons, 
can produce a relatively high yield of quasimonoener- 
getic neutrons in the forward direction via the ''Li(p,n) 
reaction. In particular, this reaction may be useful to 
provide neutrons with energies near 14 MeV, for mate- 
rials testing for the fusion program. 

In addition to the above application, quasimonoen- 



ergetic and broad neutron sources are needed for re- 
search in other areas in nuclear science and technology, 
e.g. Accelerator Transmutation of Waste (ATW), radi- 
ation damage studies, medical isotope production, and 
physics cross section experiments. 

In order to assess the feasibility of using a lithium 
target to produce neutrons, accurate evaluated cross 
section data are needed. These data can be used to 
predict the neutron yield, as well as the neutron energy 
and angular dependencies, from both thin and thick 
targets. (Thin targets allow the possibility of producing 
quasimonoenergetic neutron sources^, whereas thicker 
targets produce broad-spectrum neutron sources.) 

Table-based charged-particle transport is a new fea- 
ture of the MCNPX transport code system and has 
been implemented for protons in the developmental ve- 
sion [3]. Proton evaluations for energies up to 150 MeV 
have been completed about a year ago for 42 isotopes 
[4]. The present study extends our library for ''Li, not 
covered by [4]. A detailed report on our work may be 
found in [5]. 

11. AN OVERVIEW OF AVAILABLE ^Li(p,xn) DATA 

An overview of reactions used for the production 
of fast quasimonoenergetic neutrons, including the 
^Li(p,xn) reaction, may be found in our report [5] and in 
the comprehensive review by Drosg [6] . Natural lithium 
consists of isotopes ^Li and ^Li with abundances 7.42 
and 92.58%, respectively. The ''Li(p,n)^Be reaction was 
reviewed in 1960 by Gibbons and Newson [7]. Useful 
information on kinematics and technical aspects related 
to the use of "'^*Li may be found in [8]. 

Experimentally, the ^Li(p,n) reaction was measured 
by several groups: 

A) From threshold to about 8 MeV, by Meadows and 
Smith, using the Argonne Fast-Neutron Generator [9, 
10]; 

B) From 10 to about 20 MeV, by Anderson et at, at 



^Though our results shown in this paper indicate that once 
the incident proton energy exceeds a few-MeV, even thin targets 
result in a substantial number of lower-energy neutrons. 
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Lawrence Radiation Laboratory [11], and from 4.3 to 26 
McV, by Poppe et al. using EN tandem Van dc GraafF 
accelerator (and AVF cyclotron, for Ep > 15 MeV) at 
Livermore [12]. Poppe et al. have been summarized in 
Ref. [12] practically all measurements of this reaction 
before 1976; 

C) At 15, 20, and 30 MeV, by McNaughton et al, at 
the Croker Nuclear Laboratory [13]; 

D) By a number of other groups (see, e.g., [14-18]), 
at proton energies above 20 MeV. 

Between 1.9 and 2.4 McV bombarding energy, the 
neutrons are monoenergetic and the reaction cross sec- 
tion is large. Therefore the ^Li(p,n)''Be reaction has 
long been used as a source of neutrons (no) at these 
energies [7]. 

Above 2.4 MeV the first excited state of '^Be at 0.43 

MeV may be excited, producing a second group of neu- 
trons (ni). However, below 5 MeV the zero-degree 
yield of these low energy neutrons is less than about 
10% of the ground-state yield, so that the usefulness of 
the reaction as a monoenergetic neutron source is only 
slightly impaired. 

Above 3.68 McV, the threshold for the three-body 
breakup reaction ''Li(p,n'^He)^He, neutrons from this 
mode contribute also to the total neutron yield. The 
zero-degree neutron spectra from this mode arc very 
broad, of an evaporative-type form, and were measured 
up to Ep = 7.7 MeV by Meadows and Smith [10]. 

Above 7.06 McV, the threshold for the reaction 
^Li(p,n)''Be with excitation of the second excited state 
of ^Be, neutrons from this mode (n2) also begin to con- 
tribute to the total neutron yield, but this contribution 
is not significant [10]. Therefore, the usefulness of the 
''Li(p,n)'^Be reaction as a source of neutrons at ener- 
gies above 5-7 MeV would depend on a particular ap- 
plication. Only a very good energy resolution would 
allow one to separate the substantial number of ni neu- 
trons to obtain a monoenergetic source. However, if 
the application can tolerate including both no and ni 
neutron groups, then the reaction has the favorable fea- 
ture of use of a cheap solid target with a forward an- 
gle laboratory cross section approaching 7.4 mb/sr at 
Ep = 15 MeV and 14.5 mb/sr at Ep = 20 MeV. It 
should be noted that although the neutron yield is large, 
the presence of low energy neutrons from different reac- 
tion modes may limit the usefulness of the ^Li(p,n)''Be 
reaction as a monoenergetic neutron source at proton 
energies above 10 MeV [12]. 

The no and ni neutron yields are strongly forward 
peaked. Their measured cross section increases mono- 
tonically with increasing Ep (see Table II in [5]). This is 
not true for the total neutron cross sections integrated 
over all emission angles, shown for Ep < 26 MeV in Fig. 
1 (see [5] for a compilation of all available data). 



One can see that in contrast with the zero degree 
differential cross section, the integrated cross section 
decreases monotonically above ~ 5 MeV, reflecting the 
onset of the strong forward peaking. Above Ep — 7 
MeV, the total number of ni neutrons which leave ^Be 
in the first-excited states is always greater than 25% 
of the number of no neutrons which leave ""Be in the 
ground state, reaching a maximum of about 55% at 
Ep'^9 MeV. 

A compilation of presently available measured lab- 
oratory differential cross sections at zero deg for the 
■^Li(p,n)'^Be no and ni (0.0 MeV and 0.43 MeV levels 
of ''Be) reactions and the ratio of ni to Uq neutrons 
and of measured laboratory integrated cross sections as 
functions of proton kinetic bombarding energy may be 
found in our report [5]. 
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Figure 1: Measured total cross sections for the 
^Li(p,no)''Be and ^Li(p,ni)^Bc reactions between 4 and 
26 MeV [12] together with our evaluations in the 
ENDF-6 [19] format. 

III. EVALUATION METHODS 

The production of an evaluated data library for pro- 
ton reactions on lithium poses certain difficulties. In 
particidar, nuclear model calculations based on statis- 
tical preequilibrium and equilibrium decay theories be- 
come unreliable for light target nuclei, where the nu- 
clear levels are widely spaced. For this reason, one 
would like an evaluated data library that is based, as 
much as possible, on measured data. 

We have therefore adopted the following approach. 
There are numerous measurements of the important 
''Li(p,no) and ^Li(p,ni) direct reactions populating the 
ground and first excited states of ^Bc, and we have 
used these experimental data for the cross sections and 
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angular distributions in the evaluation. This is impor- 
tant because many aspects of the production of neutron 
sources via ^Li(p,n) reactions involve use of the high- 
energy quasimonoenergetic uq and ui neutrons. To 
model the remaining cross section available for nuclear 
reactions, i.e. the overall reaction cross section minus 
these direct reaction cross sections, we use GNASH nu- 
clear model calculations [20,21]. 

Use of the GNASH calculations is expected to lead to 
some weaknesses in the evaluated data, because of the 
aforementioned difficulties in using a statistical model 
code. However, we note that the calculations do, at 
least, include certain important constraints such as en- 
ergy, flux, angular momentum, and parity, conserva- 
tion laws. Since the GNASH calculations make use 
of experimental nuclear levels information, they can 
also be expected to lead to emission spectra that cor- 
rectly include peaks and gaps in the calculated energy- 
dependent spectra at the correct locations, i.e. peaks 
at energies where final states can be excited, and gaps 
where there are no final nuclear states available. 

A. Evaluation of no and ni cross sections and angular 
distributions 

The '^Li(p,no) and ''Li(p,ni) evaluations were based 
upon measured data. The center-of-mass measured an- 
gular distributions of both the Uq and ui neutrons were 
fitted using Legendre polynomials, 

, N 
n=0 

For proton incident energies up to 12 MeV, such Leg- 
endre fits have been published by Poppc et ai, and 
therefore we adopted their results (Fig. 5 of Ref. [12]). 

At higher proton incident energies, we performed 
Legendre polynomial fits to measured angular distribu- 
tions. For proton incident energies of 15.1, 15.9, 16.9, 
17.9, 18.9, 19.9, 20.9, 22.0, 23.0, 24.0, 25.0, and 26 MeV, 
there are available measurements separately for no and 
ui by Poppe et al. [12]. Above 26 MeV, there are 
fewer experimental data at only a limited number of 
incident energies and emission angles, and furthermore, 
only data for the sum of uq and ni are available. We 
used the measurements by Batty et al. [22] at 30 and 
50 MeV, by Goulding et al. [23] at 119.8 MeV, and by 
Watson et al. [24] at 134.2 MeV. We fitted these data 
using for the ratio R = n\/no a fixed value of 0.35, 
as measured by Poppe et al. [12] at Ep = 25.9MeV, 
and smoothly interpolated/extrapolated the resulting 
Legendre coefficients up to Ep = 150 MeV. 

As an example. Fig. 2 shows the sum of no and ni 
measured total production cross section (for references. 



see [5]) as a function of proton kinetic energy, Ep, com- 
pared with our evaluation. One can see good agreement 
in the whole energy range up to 150 MeV. 




Proton Kinetic Energy (MeV) 

Figure 2: Comparison of the sum of measured [12, 15, 
16, 24, and 25] laboratory angle integrated cross sec- 
tions for the ^Li(p,no)^Bc and ^Li(p,ni)^Be reactions 
(for details, see [5]) with our evaluation in the ENDF-6 
format. 

B. GNASH calculations 

The latest version of the GNASH code has been de- 
scribed in Ref. [20] , and its latest application in nuclear 
reaction evaluation work has been described in Ref. [4] . 
For this reason, here we provide only an overview of the 
models used in the calculations, concentrating on new 
features. 

GNASH calculations of preequilibrium and Hauser- 
Feshbach decay require input parameter information 
describing the optical model transmission coefficients, 
nuclear level densities, gamma-ray strength functions, 
as well as nuclear level information for all the nuclei 
that can be populated in the reaction. For most of these 
quantities, we used default parameter information, as 
described in Ref. [4]. Below, we provide information 
on the optical model we used for ^Li. 

Although the appropriateness of an optical potential 
for nucleon scattering on a light nucleus such as lithium 
is questionable, we have, for pragmatic reasons, made 
use of a relativistic potential for ^Li published by Chiba 
et al. [26]. This potential was shown to provide a sur- 
prisingly good representation of available elastic scat- 
tering and total cross section data, for energies up to 
a few hundred MeV [26]. We converted this neutron 
potential for use on ^Li by making small isospin cor- 
rections, and produced a version for proton scattering 
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by including Coulomb correction terms. The modified 
potentials resulted in a calculated total neutron cross 
section that accounted for measured data fairly well. 
The calculated proton reaction cross section was also 
in good agreement with measurements (see Fig. 3). 
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Figure 3: Proton reaction cross section of ^Li. The solid 
curve shows values calculated with the optical model by 
Chiba et al. [26] modified here as described in Section 
III.B. Experimental data are from Refs. [27, 28]. 



Since the GNASH code includes uq and ni neu- 
tron emission contributions from compound nucleus 
and preequilibrium decay, these cross sections had to 
be ignored since, in our evaluation, they are based on 
measurements. Therefore we wrote a utility code that 
takes a GNASH output, modifies the ng and ni neutron 
emission cross sections to zero, and in the evaluated file 
introduces the cross sections based on experiment as 
described in Sec. Ill A above. 

Examples of the ^Li(p,xn) calculated zero-degree 
double-differential cross sections of neutrons are shown 
in Fig. 4, compared with measurements. Solid 
lines show the results from the GNASH calculations; 
Dashed-lines show the no and ni contributions (evalu- 
ated from the measurements). Comparisons are made 
with Ep = 15, 20, and 30 MeV data of McNaughton et 
al. [13], Ep = 40 MeV data of Jungerman et al. [14], 
and at Ep = 55, 90, and 140 MeV data of Byrd and 
Sailor [29]. In these figures, in addition to the highest- 
energy peak due to no and ni neutrons, several lower- 
energy peaks due to n2 and na neutrons, representing 
the population of the 4.55- and 6.51-MeV states in ''Be 
are evident (see also, e.g.. Fig. 2 in Poppe et al. [12], 
where even the peaks by n4 and ns, representing the 
population of the 7.19- and 10.79-MeV states in ''Be 



are clearly observed). 

We emphasize that the GNASH results shown in 
Fig. 4 have not been broadened to account for the 
experimental detector resolution. This explains why 
peaks and dips are observed in the calculated energy- 
dependent spectra, but are not seen so prominantly 
in the experimental data. Overall, it is evident that 
the evaluation accurately represents the no and ni neu- 
trons, as it should since it is based on the experimental 
data. The lower-energy neutrons from the GNASH cal- 
culation agree with the measurements less well, but it is 
hoped that the accuracy obtained is sufficient for LEDA 
design calculations. 

IV. ENDF DATA 

In this section we show some graphical representa- 
tions of the evaluated ENDF data. These figures help 
illustrate certain features of the evaluation. The figures 
shown here are a small subset of figures produced auto- 
matically by the N JOY code [30] . The full set of figures 
can be viewed on the T-2 WWW site, by first going to 
http://t2.lanl.gov/data/hc.html, registering, and then 
clicking on prcmade plots in Postscript format for pro- 
tons (http://t2.lanl.gov/data/p-ps/). 

Figure 5 shows the proton elastic scattering angular 
distribution for incident energies up to 150 MeV. The 
increased forward-peaking with increasing incident en- 
ergy is evident. Figure 6 shows the angular distribu- 
tions, reconstructed from the ENDF Legendre coeffi- 
cients, for the no and ni emitted neutrons. 

Four separate figures are combined together in Fig- 
ure 7. On the top left, the overall cross sections for the 
no, and ni reactions are shown, as is also the "remain- 
ing" cross section, designated as MT=5, that is used 
in the GNASH calculations. The sum of all these cross 
sections is the reaction cross section as shown in Fig. 3. 

In Fig. 7, top right, the inclusive production cross 
sections are shown for the light ejectiles. Note that, 
unlike for reactions on heavier nuclei, these do not keep 
increasing with incident energy up to 150 MeV, but 
instead become approximately constant above 40 MeV. 
This energy corresponds, approximately, to the total 
binding energy of ^Li, and therefore above this energy it 
is not possible to obtain additional particle production 
by increasing the incident energy. 

In Fig. 7 lower- left, a perspective plot is shown 
of the angle-integrated neutron emission spectra ob- 
tained from GNASH, as a function of incident en- 
ergy. The increasing role of preequilibrium high-energy 
neutron emission is seen as the incident energy in- 
creases. Finally, in Fig. 7 lower-right the total heating 
(MeV/collision) is shown as a function of incident en- 
ergy. This includes energy transferred to all secondary 
charged-particles, except protons (since this is the pro- 
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jectile, and it is therefore assuraed that a transport cal- 
culation will always be explicitly tracking the secondary 
protons). 

V. SUMMARY AND FUTURE WORK 

Our new evaluation for protons on ^Li is available in 
the ENDF-6 format for use in radiation transport calcu- 
lations. By using and testing these data in simulations 
of thick and thin lithium targets, we hope to obtain 
feedback on possible improvements that are needed. 

An example of a recent test [3] of our library is shown 
in Fig. 8. We have modeled the 43-MeV proton source 
from JAERI's TIARA AVE cyclotron with MCNPX. 
Protons impinge on a 3.6-mm thick ^Li target. Re- 
sulting neutrons are constrained by an iron collimator 
10.9 cm in diameter and 225 cm long. We have mod- 
eled this target, assuming a monocncrgctic point proton 
source. The neutron flux is tallied on the surface ex- 
iting the collimator. In the report with experimental 
results [31], the neutron flux is normalized to unity in 
an energy band between 36.3 and 45.5 MeV. We have 
normalized our results in the same manner. Two MC- 
NPX calculations were performed, one using tables by 
our library and the other using the Bertini Intranuclear 
Cascade Model (INC) in conjunction with a preequilib- 
rium-|-evaporation model. 

Results for the experiment and the two calculations 
are shown in Fig. 8. Overall, neither calculation is in 
completely satisfactory agreement with the experiment. 
This is not surprising given we used a cascade, a pree- 
quilibrium, and an evaporation model in the Bertini 
INC calculations and the physics of the GNASH code 
for such a light target. However, we observe that the 
width of the neutron peak more closely matches the ex- 
periment when proton tables from our library are used. 
This could be expected in advance, since we used for 
the peak a parametrization of experimental data in our 
tables. The Bertini INC predicts a neutron peak that 
is lower in energy and much broader than observed. 

For completeness, we list here a couple of areas where 
we know there is room for improvement. The first point 
emphasizes weaknesses in the ENDF evaluation that 
may be problematic for simulations above 30 MeV for 
our case (see details in [5]). 

Secondly, our analysis was not developed to include 
an accurate treatment of gamma-ray production, but 
instead focussed on neutron emission. Indeed, the nu- 
clear level information in the GNASH calculations de- 
scribing the decay of excited states by gamma-ray emis- 
sion was not carefully checked. Therefore, if gamma-ray 
emission becomes of concern in LEDA target design 
studies, an upgrade to this evaluation may be needed. 
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Figure 4: Zero degree neutron emission from the ^Li(p,xn) reaction. Solid lines show the results from the GNASH 

calculations, dashed lined show no and ni contributions evaluated from the measurements. The GNASH-calculated 
spectra have not been broaded using the experimental resolution. Data are from Refs. [13, 14, 29]. 
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Figure 5: Proton elastic scattering angular distribution of ^Li at different incident energies up to 150 MeV. 




Figure 6: Angular distributions for uq and ni neutrons from p + ''Li interactions at different incident energies 
to 150 MeV. 
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Figure 7: Different channel cross sections, neutron angle-integrated energy spectra, and heating (MeV/coUision) 
functions of incident energy for p + ^Li interactions, as indicated. 
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Figure 8: Neutron flux from a 3.6-mm thick ^Li target normalized to unity in the energy band between 36.3 
and 45.5 MeV calculated with MCNPX using our data library (dashed histogram) and results by the Bertini INC 
in conjunction with a preequilibrium+evaporation model (long-dashed histogram) compared with the experimental 
data from [31] (solid histogram). 
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